Collisions of cold magnesium atoms in a weak laser field 
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We use quantum scattering methods to calculate the light- 
induced collisional loss of laser-cooled and trapped magne- 
sium atoms for detunings up to 30 atomic linewidths to the 
red of the 1 So- 1 Pi cooling transition. Magnesium has no hy- 
perfine structure to complicate the theoretical studies. We 
evaluate both the radiative and nonradiative mechanisms of 
trap loss. The radiative escape mechanism via allowed 1 E U ex- 
citation is dominant for more than about one atomic linewidth 
detuning. Molecular vibrational structure due to photoasso- 
ciative transitions to bound states begins to appear beyond 
about ten linewidths detuning. 

34.50.Rk, 34.10.+X, 32.80.Pj 



Light-induced collisions between cold, neutral alkali 
atoms have been widely studied experimentally and the- 
oretically in magneto-optical traps PJ-pl. Such collisions 
cause loss of atoms from the trap due to radiative or 
nonradiative molecular processes after laser excitation 
of molecular states of the atom pair at large internu- 
clear separation R. These loss processes for alkali atoms 
are still not well-understood when the laser detuning A 
is only a few natural atomic linewidths r at to the red 
of atomic resonance, primarily because of complications 
due to the molecular hyperfine structure in the alkali 
species [||3] . Molecular vibrational and rotational struc- 
ture in the spectrum of loss versus A is unresolved due to 
broadening by fast spontaneous decay and the high den- 
sity of molecular states. On the other hand, high resolu- 
tion photoassociation spectra for trap loss at much larger 
detunings are quite well-understood quantitatively, even 
for alkali species (|J^]. 

The study of trap loss collisions in alkaline earth 
species offers a significant opportunity to improve our 
understanding of collisional loss processes in the small 
detuning regime. The major isotopes of Be, Mg, Ca, Sr, 
and Ba have no nuclear spin, and thus lack hyperfine 
structure. Alkaline earth atoms can also be cooled and 
trapped (e.g. Mg §, Ca §, Sr |§). They offer an ex- 
cellent and previously rather unexplored testing ground 
for cold collision theories. A recent study of Sr has taken 
a first step in this direction ||. One intriguing aspect 
of the alkaline earth collisions is the possibility to ob- 
serve the effect of relativistic retardation corrections on 
the trap loss dynamics by exciting a molecular state with 
a dipole-forbidden transition at small R which becomes 



allowed at large R. This predicted effect for the 2 U state 
of alkali species |1(| is masked by other states, but may 
be prominent for the 1 n s state of alkaline earth species. 
In this paper we calculate collisional loss rate coeffi- 
cients of Mg atoms in a magneto-optical trap, induced by 
a weak laser field detuned up to 30T at to the red of atomic 
resonance (T at = 27T-80 MHz). This spans the range from 
small detuning, where the spectrum is unstructured, to 
the far-detuning range, where resolved photoassociative 
features due to specific molecular vibrations begin to ap- 
pear. For 24 Mg we assume that the cooling transition is 
1 So- 1 Pi, (A = 285.2 nm). Fig. |l| shows the short-range 
part of the Mg2 potentials |flj. A laser with A ~ 1- 
30 T at excites molecular states with attractive potentials 
at some Condon distance Re ^ 250 a (a = 0.0529177 
nm). In 24 Mg these states can be either the 1 Y,+ state or 
the 1 LT g state. Spin-orbit coupling -ffso mixes 1 TL g with 
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g state at a curve crossing around R 



7 at) and 
ao- After a change 
of molecular state the quasimolecular atom pair can sep- 
arate to a ^o atom and a 3 P atom ( 3 Po, 3 Pi or 3 P2), 
with a shared kinetic energy increase equal to about 0.06 
atomic units (AE/ks ~ 19 000 K). These hot atoms can 
not be recaptured by the trapping and cooling lasers and 
they are thus lost from the trap. In alkalis this process 
arises from collisional transfer between two fine structure 
states of the same electronic configuration, and is called 
the fine structure changing mechanism (FS). Here we call 
it the state changing mechanism (SC) between different 
electronic states. 

After excitation at Re the quasimolecule obtains ki- 
netic energy via internal acceleration by the attractive 
potentials. If it decays back to the ground state, it pre- 
serves this gain in kinetic energy, thereby increasing the 
kinetic energy of the asymptotic collision products (two 
4 So atoms). If the energy gain exceeds the trap depth, 
the atoms can escape from the trap. This is the radia- 
tive escape mechanism (RE). Spontaneous emission has 
a crucial role in both the SC and RE loss mechanisms, 
and has to be included in the theoretical treatment. 

The structure of the quasimolecule at large R is deter- 
mined by the dipole-dipole interaction, for which we use 
the long-range retarded potentials from Ref. |12|. These 
join with the short-range wells in Fig. (p to give poten- 
tials that support a number of bound vibrational states 
of the quasimolecule. The molecular linewidths with rel- 
ativistic retardation corrections for the two relevant ex- 
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FIG. 1. The molecular states of Mg2 in atomic units (1 
a.u. = 4.3597482xl0~ 18 J) corresponding to the asymptotic 
atomic states 1 So+ 1 Pi and 1 So+ 3 Po.i,2 (ill; the zero of en- 
ergy is at the ground state 1 So+ 1 So asymptote. There are 
four states related to each asymptote, of which two are at- 
tractive and two are repulsive at large R, where the system is 
expected to be resonant with the laser field. Each set of two 
states has one II state and one E state. 



cited states with attractive potentials are fl^] 



rs„ (R) = r at < 1 -[u cos(w) - sin(u)] 



(1) 



rn 8 (i?) = r at l 



2m 3 



[mcos(m) - (1 - u 2 ) sin(u)] 1 , (2) 



where u = R/X, and A = X/2n = 857.7 a . For small u 
we get r Eu ~ ^ (2 - u 2 /W) and Tn, =* r at u 2 /5. Note 
that the optical coupling between the quasimolccule and 
the laser field is proportional to \/Ta (A = £„,II S ). As 
detuning A increases the probability to excite the ^g 
state is reduced. On the other hand, the survival of the 
quasimolecule on the ^^Hg state is much better than for 
the 1 S U state. We expect that RE loss comes mainly from 
the ^u state, and that qualitatively the ^^Ug state and 
the 1 S U state contributions to SC loss behave differently 
as a function of A. The important question is the relative 
magnitudes of the various loss mechanisms. 

The standard small-detuning studies of trap loss ignore 
the vibrational structure of the trap loss, and consider 
the problem as a fully dynamical one, involving only the 
electronic potentials M. This is based on the expecta- 
tion that at small detunings the vibrational structure can 
not be resolved. In the time-dependent picture this sim- 
ply means that after the excitation at large R the quasi- 



molecule decays well before it completes one or more vi- 
brations. However, here the retardation effect opens the 
1 U g state for excitation but closes it quickly for decay, so 
we might expect to see vibrational structure in trap loss 
even at small detunings if the 1 H g state SC contribution 
dominates over the 1 S U state SC and RE contributions. 

The ground state potential of Mg2 is essentially fiat in 
the long-range region of excitation. Since a number of 
total angular momentum J values contribute to trap loss 
in our detuning range, we assume that effects of the un- 
known phase shifts in the ground state wavefunctions due 
to the short-range potential is removed by summing over 
J and thermal averaging. Thus, we represent the ground 
state potential by a Lennard- Jones 6-12 potential with a 
well depth of 0.0028 atomic units and an inner turning 
point of 6.23 eto UM- Since the excited ab initio poten- 
tials do not permit a calculation of vibrational levels to 
spectroscopic accuracy, the specific forms of the short- 
range excited state potentials are not important for our 
purposes of modeling the qualitative structure and mag- 
nitude of the collisional loss. We fit Lennard- Jones 3-6 
potentials to the ab initio potentials of Ref. pT| , keeping 
the long-range form fixed to its known value. Modeling 
of experimental data should ultimately permit calibra- 
tion of the unknown phases associated with the uncertain 
inner parts of the potentials. 

A complete study must include the three-dimensional 
aspects of molecular rotation. Since all collision direc- 
tions are possible relative to laser polarization, the ef- 
fect of axis rotation is to introduce the usual molecular 
branch structure of the possible transitions. The angular 
momentum J" in the ground state can only be that of 
axis rotation, or partial wave I" = J" . In the excited 
state S, 1 II molecular basis, I is not a good quantum 
number but total angular momentum J' is. The possi- 
ble transition branches are P, Q, and R, for which J 1 = 
J" — 1, J", and J" + 1 respectively The quasimolecule 
ground state can couple to the 1 S U state only by P and R 
branches, but to the 1 Ii g state by P, Q and R branches. 
The radiative coupling terms for the two electronic states 
for each branch and partial wave I" are shown in Table |. 
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TABLE I. The radiative coupling weight factors oa,b,i" 
for each electronic state, ground state partial wave I" and 
excitation branch, obtained by summing over the degenerate 
projection states. The full radiative coupling in atomic units 

is Va,b,i« = 1.16 x 10- 4 a A ,B,i»Vl (W/cm 2 )^r A (R) (a.u.). 



Since the laser field is assumed to be weak, reexcitation 
of any decayed quasimolecular population can be ignored 
and loss rates can be calculated by the complex potential 
method [ [l3|Jl4| . Additionally, each ground state partial 
wave only couples to the rotational states of the ^u 
or 1 II ff state, but these states do not couple further to 
other ground state partial waves. Thus, the weak field 
approximation permits us to truncate the problem into 
a set of independent three-state problems for each initial 
partial wave with one ground state g, one excited state 
e, and one probe state p. The total thermally averaged 
loss rate coefficient via state e is 
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B=P,Q,R 



where e is the relative asymptotic initial kinetic en- 
ergy, Qt = (27r/xfcBT//i 2 ) 3 / 2 is the translational parti- 
tion function (/x is the reduced mass of the two colliding 
atoms), and IS'gpl 2 represents the probability of a tran- 
sition to a probe channel p which simulates the effect of 
the SC or RE exit channels [Q. Identical particle ex- 
change symmetry ensures that only even partial waves 
exist for the ground state. The ground state centrifugal 
potential provides a natural cutoff £ m ax to the sum over 
Ccn m Eq. (f) where h 2 £ max (£ max + l)/2/ii? 2 > e at the 
classical turning point Rt for collision energy e. 

We calculate S gp (e,A,l,B) using a complex poten- 
tial [J13|,[l4|, i.e., spontaneous decay is taken into ac- 
count by adding an imaginary decay term — ihTA{R)/2 
(A = S„,n s ) to the excited state potentials. We solve 
the corresponding time-independent Schrodinger equa- 
tion with the appropriate boundary conditions and ob- 
tain the S'-matrix elements. 

For the SC loss studies the probe state represents 
the actual loss channel due to the curve crossings men- 
tioned above which cause transitions to the molecular 
states which separate to 1 So+ 3 P atoms. The short- 
range molecular spin-orbit couplings are approximated 
using Table Al of Ref. Q. The 1 T, U - 3 U U and ^Tlg- 3 ^ 
matrix elements are A/v2 and A/2 respectively, where 
A = 1.84 x 10~ 4 atomic units is 2/3 of the atomic 3 P2- 
3 P splitting. 

For the RE studies we assume that the trap depth is 
0.5 K, and set the probe state to cross the excited quasi- 
molecular state at the point R e where a kinetic energy 
increase of 1 K for the atom pair has been gained after ex- 
citation. Equation (5) of Ref. pi shows that \S gp \ 2 can 
be factored into two terms: IS'gpf = Jq -Pdccay, where Jq 
is the probability of being excited from g and reaching 
R e on e (Jq accounts for the possibility of multiple vi- 
brations). Jq is calculated from l-S^pl 2 as explained in 
Ref. [Q, and Pdecay, the total probability of radiative 
escape during a vibration in the region R < R e , is calcu- 
lated using a classical trajectory: 



decay 



= 1 — exp(— a), 



a = 2 f dR V J^l, (4) 
Jr b 



v(R) 



Ts„ (R) is given by Eq. (||), and v(R) is the classical speed 
related to the ^^ state potential. 

Fig. g shows calculated rate coefficients K for SC loss 
for both states at laser intensity I — 1 mW/cm 2 . K 
scales linearly with / in the weak field limit. When A 
is small, the magnitude of K for the 1 'E U state is small 
and vibrational state structure is clearly suppressed even 
before thermal averaging; some structure begins to ap- 
pear beyond A « 10r at . However, for the 1 U g state 
strongly modulated vibrational structure persists over 
the whole range of detuning with comparable peak mag- 
nitudes. The ^^Hg vibrational structure is broadened but 
not eliminated by the sum over rotational state structure 
and by thermal averaging. The specific vibrational peak 
locations versus A for both states depend on our par- 
ticular choice for the excited state potentials and may 
be expected to shift if the correct potentials were used. 
However, the magnitudes and qualitative nature of the 
oscillations will be independent of this choice. For very 
small dctunings the -"-n^ state dominates SC loss, due 
to much more favorable survival after excitation. For 
A > 10r at 1 S„ and 1 U g contributions become roughly 
equal. The survival probability relative to spontaneous 
decay at long range approaches unity for the l Ti u state 
when A > 10r at . Note that the magnitude of the rate 
coefficients for both 1 E„ and l Ii g states is rather small, 
on the order of 10 -14 cm 3 /s. This is because of the rela- 
tively small spin-orbit matrix elements at the short-range 
crossings. 

The RE loss for the 1 T, U state, shown in Fig. ||, qualita- 
tively follows that of the corresponding ^^ state SC loss 
(see the inset in Fig. ||), but is about 200 times stronger. 
Because of this inherently larger magnitude for RE loss, 
trap loss for Mg is predominantly due to the 1 H g state 
only for very small dctunings on the order of T at . Radia- 
tive escape via E u excitation is dominant at all larger 
detunings and masks any vibrational structure due to 
1 Iig SC loss. Vibrationally resolved (but not rotation- 
ally resolved) photoassociative structure should begin to 
appear due to £ u excitation for detunings larger than 
around 1 GHz, with a rate coefficient on the order of 
10~ 12 cm 3 /s for I = 1 mW/cm 2 . 

Our calculations for Mg shed some light on Sr colli- 
sions. Since the spin-orbit coupling strength A is a factor 
of 9.6 larger for Sr than for Mg, the SC loss will be much 
stronger in Sr than we find here for Mg, and thus will 
compete with RE as a loss mechanism. If we simply scale 
up our peak magnitude of 1 Tl g loss at A = 4r at linearly 
in intensity / and by the increased A 2 factor (differences 
in mass also need to be taken into account), the result- 
ing rate coefficient has the same order of magnitude as 
that which was recently measured S . This suggests that 
Sr trap loss at small detuning, unlike that for Mg, has a 
prominent contribution from SC by the l Yi g state, consis- 
tent with the understanding of Ref. M . Since it is likely 
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FIG. 2. The SC loss rate coefficient A"sc as a function of 
laser detuning A for Mg at laser intensity 1=1 mW/cm 2 
and T = 2 mK (the Doppler temperature for the chosen Mg 
cooling transition is TDo PP icr = 1.9 mK). The dotted lines 
show the sum over branches and partial waves for e — ksT, 
and the solid lines are the energy averaged result, a) 1 T, U 
state, b) 1 IT 9 state. 



that Sr will have a modulated photoassociation spectrum 
at small A analogous to Fig. Eh, measuring trap loss ver- 
sus A would be a very useful experiment to try for the 
Sr system. Such modulations might also be observed in 
both Sr and Mg by detecting, 3 Pp. 1,2 product atoms, as 
in similar alkali experiments §4|Jl6||. 

Exploring the alkaline earth systems offers the 
prospects for some very fruitful science. There is clearly 
a need for better ab initio or semiempirical potentials 
and long-range dispersion coefficients. Photoassociation 
spectroscopy over a wide range of detuning values should 
help characterize ground and excited state potentials and 
ground state scattering lengths 0] . The recent cooling of 
Sr below 1 /J.K using the 1 So- 3 Pi transition p raises the 
possibility of achieving Bose-Einstein condensation with 
an alkaline earth species. Comparing experimental trap 
loss for alkaline earth isotopes with and without hyper- 
fine structure might also shed some light on the role of 
hyperfme structure in alkali trap loss. 
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FIG. 3. The RE loss rate coefficient Kb,e for the E u state 
as a function of the red detuning A for Mg at laser intensity 
1=1 mW/cm 2 and T = 2 mK. The trap depth is set to 
0.5 K. The inset shows the temperature averaged results for 
the loss rate coefficients K for all three loss processes (short 
dash, solid, and long dash lines are from Figs. 2a, 2b, and 3 
respectively) . 
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